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Two Phases of Control System Design : 

Phase I is the design of optimal trajectories and associated inputs, that 
move a given plant from one operating condition to another, while 
minimizing some performance measure. Requires a nonlinear dynamic 
model of the specific system. 

Phase II is the design of a trajectory following controller (sometimes 
called a regulator or tracker) that provides a real-time control input 
perturbation to keep the plant operating near the designed optimal 
trajectory. Usually uses a linearized dynamic model of the specific system. 

Phase I: 

Our control philosophy is to charge the NiH2 cell in such a way that the 
damage incurred during the charging period is minimized, thus extending 
its cycle life. Requires nonlinear dynamic model of NiH2 cell and a damage 
rate model. We must do this first. 

This control philosophy is generally considered damage mitigating control 
or life-extending control. 
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Overview of NiH2 Electrochemistry 


+ 



Normal charge-discharge operation of a nickel-hydrogen cell: 


positive nickel electrode: Ni00H + H 2 0 + e 7 ^ Ni(OH ) 2 + OH 

charge 

negative hydrogen electrode: iH 2 + OH‘ H 2 0 + e' . 

charge 
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ssentialized performance model of NiH2 cell. 


storage 


diffusion 



discharge 

S + D + e ' reaction products 

charge 
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Electrode Behavior: Faraday’s Law 


For the discussion of the next two sections, consider the electrode process 


aA+bB < > cC+dD+ne 


O This equation represents both a chemical process and an electrical process. 

O The reaction rates can be completely determined using the electrical process by seeing 
that the chemical conversion can only occur if electrons are either arriving or leaving. 

O Thus, the chemical conversion rates are controlled by, or measured by, the electrical 
current passing through a given electrode. 

O Then recognizing that the rate of electron production is related to electrical current i, 
the following rate equations result : 

_ dq^_ _ c [ c id _ c dC _ a dA _ b dB 
dt n dt n dt n dt n dt 

where q e - is the charge of a single electron. 
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Electrode Behavior: Electrode Equation 


Consider the fluxes, y, of the species in the forward and reverse reactions 
at the electrode, 

i = jf - Jr. 

Assuming that the species fluxes are proportional to concentrations, yields 

• r a b j c d 

1 ~ kf C A C B ~ Kr C C C D . 

The rate constants, A;, can be related to the electrical potential across the 
electrode-electrolyte interface using free energy considerations, 

kf = koe (F'«T )( l-a )( v-v { o) 

£ _ fc e (F/RT)(-a)(v-v 0 ) 

Inserting these gives the electrode equation, 

i - F Lab (F/RT) (\—cc) (v— v 0 ) _ c d (F I RT)(-a){v-v Q )\ 
l — At 0 \C a L b C l C L D e / « 


The approach used is often referred to as the Butler-Volmer approach. 
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Linear-in-the-parameters Electrode Equation 


e propose a linear-in-the-parameters approximate solution to the electrode equation; 



where the &’s are parameters to be determined from data. 

This represents a compromise between the Tafel and the Nernst solutions of 
the electrode equation: 

The Nernst solution assumes that the current is so small as to be negligible 


v = v 0 + 


RT 

nF 


f „ A 


In 


\ C d J 


The Tafel solution assumes that the current is large in one direction or 
the other, which means that one of the two exponential terms is negligible 


v = v 0 + 


RT 

anF 


RT 


ln(*o) -ln(/) 


anF 
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Essentialized Model Overview 


Terminal behavior: 


current into the battery is +i, 

the terminal voltage is +v, 


stored material with self-discharge: 


*- <0 ±c, m 


dt 


R 


sd 


diffusing material: 


C d ( 0 C bulk ^d A q i{t) ? 


electrode-equation: 

v (c s , Q , 0 = k l +k 2 ln(l + z l)sgn(z') + k 3 In (c d ) + k 4 ln(l - c, ) 
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Parameter Determination 


The cell chosen is the NSWC Crane Pack ID 3602G (Gates): 
rated at 65 AHr 
uses 3 1 % KOH concentration 
maintained at 1 0 degrees C 
charge-discharge profile is a square wave 
with 35% depth-of-discharge (DOD) 

1 04% recharge ratio 

current is 26.29 A for 54 minutes charging 
-37.92 A for 36 minutes during discharge 
note that 65 AHr = 3900 AMin, 35% of 65 AHr = 1365.1 
AMin. 
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Essentialized Model with Identified Parameters 


Terminal behavior: 


current into the battery is +/, 

the terminal voltage is +v, 


stored material with self-discharge: 


dc s (t) 

dt 


= *•(0-0.0002085^(0 


diffusing material: 


c d (t) =1 - 0.00 1036 ,/// jC)034 ;'(f ) 


electrode-equation: 


v 


= 1 .3656 + 0.0265 ln(l + 1 i |)sgn(0 + 0.0229 ln(c^ ) 

- 0.0262 ln((l. 005 * 3900 - )/3900) 
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voltage 


voltage vs time for least squares identified model, cycle 10000, from akronbatt Is.m 
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A typical charge-discharge cycle, charging. 
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Damage Mechanisms for NiH2 Batteries 

Formation of y- phase NiOOH: 

overcharge 


P -NiOOH 


> y -NiOOH 


charge 


A 

A 

i 


discharge charge 


1 

r 

T 


discharge 


P -Ni(OH): 


a -Ni(OH): 


standing in KOH 

Bode’s Solid phase relationships for a NiOOH electrode. 
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Damage Mechanisms for NiH2 Batteries 


Formation of 02: 

Overcharge: 1) continuing to charge the cell after all the p -Ni(0H)2 

has been converted 
2) the charging current is too large 

The effect of this is the formation of O 2 at the nickel electrode, 
along with heating. 
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Damage Mechanisms for NiH2 Batteries 


Heating: 

Sources: 1) heat of reaction 

2) formation of O 2 

3) electrical current 

Results in: 1) the formation of r-NiOOH: 

which a) reduces cell capacity 

b) does physical damage to the cell 

Results in an increase in self-discharge reaction rates. 
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Battery Continuum Damage Modeling 


• Many Possible Damage Mechanisms 
•Hard to model all these 


• Overall Birth to Death Data will be Used Instead 


• Crane Database Provides much Information 


• Green-Hoffman Data Taken as Starting Point 
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Battery Continuum Damage Modeling 


Green-Hoffman DataNiH 2 atT = 10 °C 


1 

DOD 

Cycles 

To 

Failure 

35% 

38,000 

50% 

19,000 
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Green-Hoffman Battery Life Model 


* io 4 



DOD (%) - Depth of Discharge 
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Continuum Damage Model 


Based on G-H Data 


N f =1 885.04e 4 ' 621(1_DOD) 

JGH 


= 19151 \.13e~ 4 - 62WOD , at 10° C 
For constant damage per cycle 

D = — — = 5.222 xl0“ 6 e 462lDOD 

^ Nr 

JGH 

DOD =c x v a 

D = 5.222x1 CT 6 e 4 - 621c ' Va 

eye 
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Continuum Damage Model 


f S(v)dv = D cyc = 5.222 xlO -6 e 4 ' 621ClVa 

cycle 


where 



voltage referred 

damage rate 


For damage on charging only 

v 

*max 

f 8{v)dv = D cyc = 5.222x 1(TV- 621ciV « 

V 

*min 
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Continuum Damage Model 


v — v * 

^max y mm 


f S(v + v min )dv = 5.222xl0~ 6 e 4 ' 621ClV “ 


o 


= 5.222x10“° e 


^ f’ max ^min ) 


Thus inferring 


S(v + v : ) s 5.222 x 1 0 -6 4.62 1 c,e 


4.62 lq v 


= 2.4131x10 q e 


-5 _ _4.621qv 


Hence 
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Continuum Damage Model 


Instantaneous Damage Rate 

• . . dD dD dv cv / w • s \ 
D(t ) = — = —— = S(v(t))v(t) 

dt dv dt 


Requiring Positive Damage 


m = 


dD dD 


dt 


dv 


dv 


dt 


= s(v{t)) kol 


D{t) = 2.4131x10 q e 


-5„ „ 4 - 621 C 1 ( v_v min) 


v(t) 


D(t) = 2.4131x10 _5 c, e 


4.621 q (v— 1.2) 


v(t) 
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Modified Continuum Damage Model 

For zero damage at zero POD 

D C yc = c 2 (e C3D0D -e C}0 ) 

D =cAe C}D0D -l)= — 

eye 2 \ ) N 

c 2 and c 3 are determined to match G-H Data 

D = 1.0404 xlO" 5 ( e 3 - 602DOD -l)=— 

N f 

Repeating the previous process gives 
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Battery Life Models 


x 10 4 
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Modified Continuum Damage Model 


In terms of current 



clD 

dv 


dv 

dq 

= S(v(tj) 

dv | 

dq 

dt 


dq 



where _ is slope of charging curve 

dq 

For zero damage when voltage rate goes negative 

D(t) = 3.7475X 10 -4 c, e X602ci (v_L2) v(t), v(t) > 0 
D(t) = 0, v(t) < 0 
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Two Phases of Control System Design : 


Phase I is the design of optimal trajectories and associated inputs, that 
move a given plant from one operating condition to another, while 
minimizing some performance measure. Requires a nonlinear dynamic 
model of the specific system. 


Phase II is the design of a trajectory following controller (sometimes 
called a regulator or tracker) that provides a real-time control input 
perturbation to keep the plant operating near the designed optimal 
trajectory. Usually uses a linearized dynamic model of the specific system. 


Phase I: 

Our control philosophy is to charge the NiH2 cell in such a way that the 
damage incurred during the charging period is minimized, thus extending 
its cycle life. Requires nonlinear dynamic model of NiH2 cell and a damage 
rate model. Now that we have this we can begin the control design process. 
The specific control philosophy is employed is generally considered damage 
mitigating control or life-extending control. 
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Performance Measure 


The performance measure to be minimized is the accumulated 
damage per recharge cycle: 

J = W fs{ C *( t f)- C s( t f)J + j^P- dt 

o at 

* , X 

c s Vf ) is the desired stored charge at the end of charge, 

w fs is the cost weighting , 
dD(t) 

—ft is obtained from the damage model, 
tf for our problem is 54 Min. 
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65%-100% recharging, voltage profile 



Comparison of 65%-100% charging methods, voltage 
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65%-100% recharging, damage profile 


x IQ ' 5 



Comparison of 65%-100% charging methods, damage. 
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65%-100% recharging, current profile 



Comparison of 65%-100% charging methods, current. 
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recharging, voltage profile 



Comparison of 60%-95% and standard charging methods, 

voltage. 
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60%-95% recharging, damage profile 


x 10' 5 



Comparison of 60%-95% and standard charging methods, 

damage. 
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60%-95% recharging, current profile 



Comparison of 60%-95% and standard charging methods, 

current. 
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Optimal Charging Summary 



Damage 
per Cycle 

Cycles to 
Failure 

% Life Extension 
relative to constant 
current 

% Life Extension 
relative to constant 
current-plus-taper 

Constant-Current 

Charging 

0.000039269 

25465 

0% 

not applicable 

Constant + Taper 
Charging 

0.000026316 

38000 

49.22% 

0% 

Life Extending 
Charging 
65%-100% Cycle 
of Figure 5.4 with 
abs(dv/dt) damage 
rate 

0.000020780 

48123 

88.98% 

26.64% 

Life Extending 
Charging 
60%-95% Cycle 
of Figure 5.5 with 
abs(dv/dt) damage 
rate 

0.000019535 

51190 

101.02% 

34.71% 

Life Extending 
Charging 
65%-100% Cycle 
of Figure 5.6 with 
only +dv/dt damage 
rate 

0.000022080 

45290 

77.85% 

19.18% 


Comparison of damage for various charging methods, assumes damage only during charge. 
Percentage life extension is reduced proportionately for damage during discharge. 
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Control Design - Phase II: Tracking 


Input Measured output 



Real-time observer structure. 
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Control Design - Phase II: Tracking 



Advanced control system using advanced real-time observer. 
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Summary 


• Control Philosophy 

• Essentialized Model Development 

• Damage Model 

• Optimal Life-Extending Charging 

• Tracking Controller 

• Real-time Parameter Identification Development 

• Application to Lithium based cells 
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